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Abstract 

Experimental analyses of charged current deep inelastic charm production - as ob¬ 
served through dimuon events in neutrino-iron scattering - measure the strangeness 
component of the nucleon sea. A complete analysis requires a Monte Carlo simulation 
to account for experimental detector acceptance effects; therefore, a fully differential 
theoretical calculation is necessary to provide complete kinematic information. We 
investigate the theoretical issues involved in calculating these differential distributions 
at Next-Leading-Order (NLO). Numerical results are presented for typical fixed target 
kinematics. We present a corresponding FORTRAN code suitable for experimental NLO 
analysis. 

1 Introduction 

Recent sets of global parton distribution functions (PDFs) Jlj-Q have reached a sufficiently 
high level of accuracy that quantization and propagation of statistical errors have become 
important issues |^. It is, therefore, even more unsettling that the strange quark PDF, 
s(x, Q^), remains a mystery P] without a fully consistent picture emerging from the compar¬ 
ative analysis between neutrino and muon structure functions 0-^, opposite sign dimuon 
production in i/Fe-DIS |]8|-[l6|, or the recently measured parity violating structure function 
AxFs 00- Given the high precision of the non-strange PDF components, this situation 
for s{x, Q"^) is unacceptable both in terms of our understanding of the nucleon structure, and 
for our ability to use precise flavor information to make predictions for present and future 
experiments. 

For extracting the strange quark PDF, the dimuon data provide the most direct determi¬ 
nation. The basic channel is the weak charged current process us —> fi~cX with a subsequent 





charm decay c — fi^X'. These events provide a direct probe of the shh-vertex, and hence 
the strange quark PDF.[| For this reason, hxed-target neutrino dimuon production will pro¬ 
vide a unique perspective on the strange quark distribution of the nucleon in the foreseeable 
future. Besides, HERA provides a large dynamic range in for the CP conjugated process 
e~s —*• VC which is valuable for testing the underlying QCD evolution Within the 

HERMES experimental program |2^ the flavor structure of the polarized and unpolarized 
sea are studied from semi-inclusive DIS where DIS Kaon production has obvious potential 
to probe strangeness. Thus, HERA and HERMES can complement fixed-target neutrino 
dimuon data with information at different energies and from different processes; therefore, 
neutrino DIS serves for now as an important benchmark process to perform rigorous and 
rehned comparisons between the experimental data and the theoretical calculations. In the 
long run, a high luminosity neutrino factory could, of course, considerably raise the accuracy 
of present day information from v DIS H- 


The theoretical calculations of inclusive charged current charm production have been 
carefully studied in the literature |p|, p!8|, p2|-p7[| . Additionally, the charm fragmentation spec¬ 
trum has also been calculated in detail ||l6|,p!8|. While inclusive calculations are sufficient 
for many tasks, a comprehensive analysis of the experimental data at NLO requires addi¬ 
tional information from the theoretical side. In charged current v-Fe charm production, the 
detector acceptance depends on the full range of kinematic variables: {x,Q^, z,ri}.^ Here, 
X is the Bjorken-x, Q is the virtuality of the hF-boson, 2 ; is the scaled energy of the charm 
after fragmentation, and rj is charm rapidity. The theoretical task is, mutatis mutandis, 
not too different from the extraction of the neutral current charm structure function F 2 as 
performed by the HERA experiments this analysis uses the theoretical calculation of 


the differential cross section ||^ to extrapolate into regions of poor acceptance. 


The organization of this paper is as follows. In Section ||, we discuss the key factors that 
influence the acceptance of the experimental detector. In Section we review the theoretical 
calculation of the fully differential cross section at NLO in QCD [|T^]. In Section ^ we present 
numerical results for typical fixed target kinematics. Finally, in Section |], we draw our 
conclusions. 


^In contrast, single muon production only provides indirect information about s{x,Q^) which must then 
be extracted from a linear combination of structure functions in the context of the QCD parton model. 

^Azimuthal ^-dependence is controlled to be flat by rotational symmetry around the axis of the boson- 
target frame. It is, therefore, not indicated. 
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2 Experimental Environment: uFe DIS 


Dimuon events from neutrino charged current charm production can provide a clear set of 
events from which to study the strange sea. Their signature in a detector is a pair of oppo¬ 
sitely charged muons and a hadronic shower originating from the same vertex. The second 
muon is produced in the semileptonic decay of the charmed particle. Detector geometry 
and applied cuts affect how these events are reconstructed, and must be corrected for when 
making measurements of the underlying charm production. These corrections require a cross 
section differential in variables in addition to the typical E,y,x, and or y of the charged 
current cross section. 


To seperate dimuon events from backgrounds a minimum energy requirement must be 
applied to the muon from charm decay. For the muon to be visible at all, it must hrst be 
energetic enough that it travels further in the detector than the particles which make up the 
hadronic shower. The background from muons from nonprompt decays of pions and kaons 
within the shower is large at low energy, and must also be reduced. Typically a cut on the 
charm decay muon’s energy of a few GeV is applied to guarantee it be reconstructable, and to 
reduce the nonprompt backgrounds. The energy of the decay muon depends on the energy 
of the charmed meson from which it decays, and therefore depends on the fragmentation 
parameter z. In NLO, it also depends on the rapidity of the charmed quark, which maps onto 
the IF-parton CMS scattering angle. A charmed meson with low ; 2 , and/or low rapidity will 
be less likely to pass a cut on the decay muon’s energy than one with high z and maximah^ 
rapidity. In addition, the linear construction of typical hxed target neutrino detectors can 
adversely affect the reconstruction of decay muons with large scattering angles. Considering 
these effects, it is important that the charm production cross section’s dependence on both 
and charm rapidity be understood to be able to model dimuon events and their detector 
acceptance properly. 


Fig. I illustrates the effect of a decay muon energy and angle cut on ^ and rapidity 
acceptance at a particular Ei^,x, and point. A toy monte carlo was constructed that 
generated charmed mesons with flat distributed z and charm rapidity, and decayed them into 
muons following |^. Detector smearing effects were not included. The hgure shows muon 
acceptance after a 5 GeV cut on energy, and a maximum scattering angle cut of 0.250 radians, 
as a function of 2 ; and charm rapidity for E^, = 80 GeV, x = 0.1, and = 10 GeV^. The 
acceptance is a smooth function in both kinematic variables, but flattens out as ^ approaches 
1 and rapidity approaches its maximum value. To properly correct for acceptance the cross 
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Figure 1: Relative acceptance A{z,ri) induced from typical kinematic cuts (as detailed in 
the text) on the decay-muon for = 80 GeV, x = 0.1, = 10 GeV^. 

section must be differential where the slope is nonzero, but in regions where the acceptance 
is flat, the dependence on these variables may be safely integrated out; i.e. we conclude that 
the integrated acceptance correction 

Jdzdn ( 1 ) 

does not resolve any local deviations of da around some point ZQ,r]Q within a range where 
[dA{z, r])/dz\ x{z — Zo) and [dA{z, h)/'^h] ^ {v~Vo) are smaller than the typical MG accuracy. 
In our example, these conditions are safely met within a range of, say, 1 unit in rapidity and 
0.2 units in away from the phase space boundary. 

3 Differential Distributions at NLO 


Recorded charged-current charm production rates must be corrected for the detector accep¬ 
tance which, as discussed above, depends on the full range of kinematic variables: {x, Q^, z, 77 }.^ 
Therefore, we must obtain the NLO theoretical cross section which is completely differential 
in all these variables. 

While X and are fixed by the energy and angle of the scattered muon, the variables 
and f] relate to the center-of-mass-system (GMS) phase space of the hard partonic scattering 
event, which at Leading-Order (LO) is W^s' c. At NLO, the quark-initiated process 
(W^s' c) receives virtual and real corrections, and we encounter a new gluon-initiated 
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process, W^g —>• cs'. Neglecting bottom contributions, Q we denote the CKM-rotated weak 
eigenstate with s', which is dehned as: 


s' = IK, cl' S+|K,e|' d 


( 2 ) 


Similarly, we will use 

K = ^7(|K,c|' + |K,/) (3) 

to denote its QCD evolution partner, he., ds'/dluQ^ = s' ® Pgq + g' ^ Pqg. 


The NLO 2-particle partonic phase space in iV = 4 -|- 2e dimensions (appropriate for a 
single charm quark of mass rUc, plus a massless gluon or strange quark) is: 


dPSo = 


1 s 
Svr 


m: 


r(l -|- e) [ dvrs 


m 


2\2 


[d(l - y)T dy 


( 4 ) 


Note that there is only one independent variable;^ In Eq. (j^), this has been chosen to be 
y = {1 + cos9*)/2, with 9* being the IK^-parton CMS scattering angle. 


The expression for the phase space in Eq. (|^) leaves us two options for computing the 
cross section: i) we can either integrate over the entire phase space, or ii) we can compute a 
singly-differential distribution. 


3.1 Singly Differential Distributions 


We can construct any singly-differential cross section as follows. If the variable S stands for 
any kinematic variable that can be expressed through y and partonic CMS energy s, we can 
easily obtain the differential cross section with respect to H with the appropriate Jacobian 
via the relation: da/d'E = (da/dy) {dy/dE). As above, y = {1 + cos6**)/2 and the partonic 
CMS energy is given by: 


s 


91 


i-r + ^ 


( 5 ) 


with dehned below in (P). In our case, we consider a distribution singly-differential in (, 
the partonic fragmentation variable 


^ Pc-Pn ^ y{s - ml) + mg 
~ q- Pn s 

^We can safely neglect the bottom quark-initiated contributions as these are suppressed both due to the 
bottom-quark mass, and also by the Cabibbo-Kobayashi-Maskawa (CKM) matrix element. 

^To count the degrees-of-freedom, we have 2x4 momentum coordinates, minus 2 mass shell conditions, 
minus 4 energy momentum conserving conditions, minus 1 rotational symmetry around the CMS axis. 
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which reduces to the scaled energy of the unobservable (colored) charm quark, ( = Ec/Ew, 
in the target rest frame. After fragmentation the observable (colorless) charm hadron (col¬ 
lectively labeled a D-meson) carries some smaller energy = E^/Ew < (. 

The unobservable variable ( is convoluted with a fragmentation function to compute the 
observable CC DIS hadronic charm structure functions: 


D^{z) 

/§ ^ ./max(2,Cmin) ^ 


+ 


27r 












(7) 


In Eq. (1^), Dc{z) is the fragmentation function (often parameterized in the form of Peterson 
et al. [^), is the scaled charm hadron energy, ^ = x(l -|- m‘l/Q'^), Cmin = ml/s = 

(1 - -^)^7(1 - -^7) and A = QV(g 2 


The outer convolution integral over the variable 

‘^Ps,g - q \ Q^) S + 


( 8 ) 


folds in the quark and gluon PDFs. We take the factorization scale [ip used in the PDFs 
equal to the renormalization scale fip used in asipn). We find it convenient to normalize 
our jFf relative to the conventional CC DIS structure functions E^ as follows:jFj^ = E^, 
JF| = F|/2, JF| = ^ 2 / 2 ^. We adopt the standard normalization of semi-inclusive structure 
functions as implicitly defined through 


Fi{x,Q^) 



E,{x,Q^,{E}) 


(9) 


with {S} any set of kinematic variables. 


The form of the hard-scattering coefficients, t)e found in |1l6|,p!8|. These hard- 


scattering coefficients, are mathematical distributions containing both Dirac-(5 and 

(generalized) plus distributions arising from the soft singularities and mass singularities 
present in the NLO QCD calculation.^ To map these distributions into C-numbers, smooth 
test functions are required; normally, it is the parton distribution functions and the fragmen¬ 
tation functions that serve as these smooth test functions. Convolutions over both and ( 

^In this article, the term distributions may denote either kinematic or mathematical ones; the individual 
meaning will be clear from the context. 
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are necessary to guarantee C-number results for the hadronic structure function on the left- 
hand-side (LHS) of Eq. O; the convolution over folds in the PDFs, and the convolution 
over C folds in the fragmentation function. 


In light of the above discussion, we clearly see that if we require structure functions 
which are fully differential in all four variables, such as z,r]) which is required to 

compute the experimental detector acceptance, we will be left with unphysical Dirac-5 and 
plus distributions. Let us note that we are facing a new problem that is intimately related 
to the single-charm kinematics of charged-current DIS. There is no direct lesson to be learnt 
from the pair-production kinematics of the well-understood neutral-current case where 
the phase space has more degrees of freedom. We now turn to address this question. 


3.2 Fully Differential Distributions 

Suppose we wish to compute the structure function which is differential in four variables. As 
a specihc example, we will choose where p^* is the transverse momentum 

of the charm particle in the boson-parton CMS frame. From Eq. (^) and Eq. 0 we will 
hnd contains mathematical distributions of the form: ~ Pt*) 

1/(Pt* - Pr*)+> where = ^{1 - C)(Cs — m?). The 5-function distributions arise from 
the diagrams that have LO (2 —> 1) kinematics, whereas the plus-distributions arise from 
diagrams with NLO (2^2). 

In principle, C-numbers are ultimately obtained from the jFf(x, distribution 

after convolution with a smooth detector-resolution function in a MC program. While such 
an approach may work in principle, there are many inherent difficulties trying to interface 
mathematical 5-function distributions and plus-distributions with a complex detector simu¬ 
lation MC program. 

In the next sections, we hrst look at the source of the singularities which appear in the 
theoretical calculation, and then hnd a means to extract the relevant distributions in an 
expeditious manner. 


3.3 Resumming Soft Gluons and the Sudakov Form Factor 

As mentioned above, if we compute the fully differential jFf (x, -s, p), or equivalently 

the hnal state is over-constrained and this yields an unphysical delta- 
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function distribution in the observable cross section. While such a distribution will yield sen¬ 
sible results for any integrated observable, the full differential distribution da/dxdQ^dzdpr 
will contain unphysical singularities. 

This affliction is the limitation of our fixed-order perturbation theory; were we able 
to perform an all-orders calculation, such singular distributions would never occur—as is 
appropriate for a physical process. In our particular case, the singular behavior of our fixed- 
order perturbation theory would be resolved were we to include soft gluon emissions as in 
the familiar Sudakov process. We will briefly review the Sudakov resummation of soft gluon 
emissions off massless particles, and then assess how best to work with z,ri) for 

non-zero charm mass. 


One of the early resummations of the QCD logarithms of the form InQ'^/q^ arising from 
soft gluon emissions was performed by Dokshitzer, Diakonov and Troian (DDT) |^. They 
obtained the leading-log Sudakov form factor which can be written generically as: 


S{qT,Q) = exp 



as{p) 
P? TT 



where qt = —prjz is the transverse momentum of a massless particle produced in DIS 
divided by fragment at ion-z. The net result of this Sudakov form factor is to smear the pt 
distribution from a singular delta-function, to a more physical form-often taken to be a 
Gaussian.^ 


The phenomenological approach we will take in this paper is divided into two steps. 


1. We regularize the NLO calculation of J^f{x, Q^, z, rf) to provide a numerical distribution 
free of delta-functions and plus-distributions. 

2. This result is input to a Monte Carlo (MC) calculation where the additional effects, 
including iterated soft gluon emmisions, is modeled by a Gaussian distribution that 
has been fit to data. 

®There are a number of different implementations of this formalism. In the Collins-Soper-Sterman (CSS) 
formalism [Q, this process is separated into a perturbative and a non-perturbative Sudakov form factor. The 
perturbative Sudakov form factor represents a formal resummation of the soft-gluon emissions. The 

non-perturbative Sudakov form factor [S^^) is introduced to deal with the singularities in the perturbative 
coupling at at /i = Aqcd 0il- Separately, Contopanagos and Sterman |Q developed a Principal 

Value Resummation which provides a prescription for regularizing this singularity. More recently, Ellis and 
Veseli derived a technique that operates purely in qt without resorting to a Bessel transform to b-space. 





To ensure that we have an accurate representation of the differential distribution suitable 
for evaluating the experimental acceptance and extracting the strange-quark PDF, we need 
to verify that the net smearing due to step 1) and step 2) combine to yield the smearing 
measured by experiment. 


The transverse momentum of charmed particles about the direction of the charm quark 


has been measured by the LEBC EHS collaboration and by the neutrino emulsion 
experiment FNAL-E531 They £t to a Gaussian distribution: 

dN 


dpq 


^-PtKPt) 


with {p"^) on the order of ~ 900 MeV. 


This observation implies that the details of the specihc regularization procedure are 
inconsequential; they can simply be compensated by adjusting the smearing of the Monte 
Carlo so that the net smearing (regularization plus Monte Carlo) yield the experimental 
transverse momentum distribution. Effectively, this parameterizes the gluon resummation 
of the Sudakov form factor. In practice, we can proceed along items a) or b) of Section 
|3^ and adjust the regularization parameters by either varying the (p^) in the G function 
in Eq. ([loD below, or by varying the bin width differential distribution. This cross-check 
ensures that we have an accurate representation of the differential distribution suitable for 
evaluating the experimental acceptance and extracting the strange-quark PDF. 


We now discuss the implementation of this procedure in the following section. 


3.4 Implementing the Fully Differential Distribution 

As demonstrated above, the singular distributions are nothing but an unphysical artifact of 
regularized perturbation theory; for any physically observable quantity, these singularities 
will be smeared by soft gluon emission to yield physical C-number distributions As 

the theoretical machinery of soft gluon resummation is not fully developed for semi-inclusive 
DIS heavy quark production, we will use a phenomenological approach.[] There are two 
viable options that we have investigated. 

a) Gaussian pT^-smearing: We may smear the singular peak which appear in the dif¬ 
ferential structure functions. 

^Note, the application of the Sudakov resummation to the case of semi-inclusive DIS heavy quark pro¬ 
duction is under current investigation; such an approach would be of use in a future very-high-statistics 
experiment where the detailed pT-distribution could be measured. 
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b) Binning: We may map the singular distributions onto C-numbers distributions by 
integrating over bins which reflect the flnite resolution of the experimental detector. 

We examine these possibilities in order. 


3.4.1 Gaussian pr-smearing 


We can implement Gaussian p^-smearing by replacing the delta function h(py* — Pt*) with 
a narrow Gaussian distribution 5p^*) centered about the perturbative value p^^*, 

and with a width of 5py*:3 


-2(P^ - 


/^/ C,* C,* r C,*\ 

G[pi ) = 


c ,*\2 
Pr^ ) 


6p: 


c,* 

'T 


TT 


( 10 ) 


Once we have a smooth distribution in p^*, we can trivially switch to any other variable, S, 
by applying the proper Jacobian:^ 


G(S,S,hS) = G{pr,Pr.Spr) X ^ (11) 

For non-zero Sp^*, we pass G-numbers to the MG. In the limit Sp^* 0, we have G{p'^*,p'^*, 0) 
^{Pt* ~ Pt*)^ original mathematical (singular) distributions are recovered. In prac¬ 

tice, a Gaussian smearing reagularization of the in Eq. corresponds to a replacement 


Wf“K'.C.yF.A) ^ Hf A, >■) = X G(p‘t'.Pt‘) 


( 12 ) 


with p^* = p^*{i' 1 C) understood as a function of and C. 


3.4.2 Finite Detector Resolution and Binning Distributions 

An alternate approach is to map the singular distributions onto G-numbers distributions 
by integrating over bins which reflect a flnite resolution of the detector. Additionally, this 
technique can avoid passing negative weights to the MG if the bin size is sufficiently large 
enough to allow the KLN theorem to effect its cancellations. This point will be further 
discussed in Section 

®Note, the actual normalization of the Gaussian is more subtle than shown here as the integration in 
will be constrained by kinematic limitations. This complication is one reason that we do not implement this 
approach in practice. 

®Note, G(S, S, (55) will not necessarily be a Gaussian. 
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In practice, the binning means that in Eq. (|^) we make the replacement 


Hf’Oii', A) ^ c, A A) = c, A) X OhIJ', C) 


( 13 ) 


where the step function G>e(Cj() = 1 if S = is inside the bin AH over which to 

integrate, and zero otherwise. Note, this means that 0s(^^C) is subject to the distribution 
prescriptions of the original as given in Refs. ||^|T8|; he., 0s(^',C) must be treated 
with care as a multiplicative factor to the non-singular “test-function” part of the 


3.4.3 Practical Considerations 


We will briefly comment why we favor implementing the binning of choice b) as opposed to 
the Gaussian pT’-smearing. For the present case of the charm transverse momentum, H = 

(or equivalently, the charm rapidity H = p), Gaussian p-p-smearing leads to a integrable but 
delicately singular behavior near the boundary of phase space - which corresponds to the 
singular point of the distributions - unless we smear the p^* variable with a width 5py* as 
large as ~ rUc. 


In the context of resummation of soft gluon emissions into the Sudakov form factor, the 
width 5p^* should be tied to a scale characteristic of the soft processes within the hadron, 
possibly with a logarithmic growth [^. For hxed-target neutrino DIS, p^* ~ rric represents 
a comparatively large transverse momentum. Additionally, there is the complication of the 
normalization of the Gaussian given the constraints of the phase space boundaries. 


Gonsequently, while the Gaussian p^-smearing may have a simple intuitive interpretation, 
we hnd the second method of using hnite width bins to be more economical in terms of GPU 
time and numerical stability. Results for this approach with H identihed with charm rapidity 
r] will be presented in the next section. 


4 Numerical Results 

We have previously addressed in Section ^ the complication of mapping mathematical dis¬ 
tributions onto G-number functions by the introduction of bins. We also encounter large 
and negative Sudakov logarithms close to the phase space boundary where, at hxed-order, 
soft single-gluon emission is enhanced. Let us, for illustration, consider the convolution of 
a plus-distribution 1/(1 — x)+, and a sufficiently smooth and monotonically decreasing test 
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function f{x) > 0. Then, the integral 

r * /w = r,, Mxjw ^ ^(1, ( 14 ) 

will be positive for Xmin not too close to 1. As in this simplified “toy” case, the Sudakov 
logarithms near the phase space boundary diverge in the limit of zero bin-width [xmm —^ 1 
in (Eq. (|TD)]; as we increase our resolution via narrow binning, we begin to resolve the 
unphysical ^-functions and “plus-distributions.” Conversely, by using broad bins, we are 
effectively integrating over enough phase-space so that the KLN theorem ensures that we 
obtain positive physical results. 


4.1 Kinematics for Rapidity Distributions 


We now compute the normalized differential charm production cross section: 

-1 -1 


da 


da 


{x,y,z,r]} — 


dx dy dz drj 


2GIMmE, 


TT (l + 


Q2 


(15) 


We can compute dcTja, ^ ^from the master equation for the differential structure function, 
Eq. (^, with the binning procedure dehned by Eq. (|^. For our variables, we choose the 
set {x,Q,z,ri} where rj is the charm rapidity evaluated in the collinear {p±^w = 0) target 
rest frame. 


In the partonic center of mass system, the charm rapidity is dehned by: 

* 1 , e:+pI^ 

p = - In ■ 


and the other necessary variables are: 


e: = 


2 E*-pl^ 


s + 
2v^ 


(16) 


Pl,c = 


s — m" 


cos 6* 


cos 9* = 1 — 2 


2V) 

(1-0(1-A?') 


(17) 

(18) 
(19) 


1-0 

We hnd the charm rapidity to be a convenient variable as the relation between the rapidity 
in the partonic center of mass system, p* and the rapidity in collinear target rest frame p 
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Figure 2: Binned differential distribution for CC neutrino-production of charm on an 
isoscalar target; the kinematics shown are for a typical wide-band beam on a hxed target: 
= 80 GeV, a: = 0.1, = iq GeV^. 


are related by a simple relation: 


V 


ri* -\— In 
' 2 






X 


— X 


( 20 ) 


We will use the collinear target rest frame as our reference frame since it is related to the 
laboratory frame by a simple spatial rotation which depends on an event-by-event basis on 
the kinematics of the leptonic vertex. 


4.2 Tracking Down Sudakov Logarithms 

In Fig. (H) and Fig. (j^) we present results for kinematics typical of a wide-band neutrino 
beam on a hxed target: Ey = 80 GeV, x = 0.1, = 10 GeV^. We plot da in 2-dimensions 

vs. z and r] in Fig. (^, and in Fig. we show 1-dimensional slices for distinct bins in ^ 
plotted against rj. Gomparing Fig. (H)-a and Fig. @-b, we see the effect of the narrow- and 
wide-binning. 
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Finally, in Fig. (H), we show results for das^x,y,z,r),Q'^} as defined by Eq. ( p!5D for the case 
where either rj or z is integrated out; again, we display this for hne-binnings of 1x100 or 
100x1, and broad-binnings of 1x5 or 10xl.[^ 

As expected, the negative Sudakov logarithms occur at large 2 ; where we have rj < 
in the integrand]^ This behavior is clearly evident in Fig. (j^) where we see that in the case 
of hne binning (Fig. (||)-a) and large (dotted curve, corresponding to 2;e[0.95; 0.96]), the re¬ 
distribution turns negative in the lowest rapidity bin. In contrast, in the case of broad binning 
(Fig. (|)-b) and large z (dotted curve, corresponding to ze[0.8; 1]), the re-distribution remains 
positive throughout the full 2 ;-range. Similarly, the ; 2 -distribution in (Fig. (^-a) is rendered 
positively dehnite if the binning is broadened at large-^. Hence, we observe that negative 
weights can be easily avoided by using sufficiently broad bins in a reasonable broadening of 
the binning 

To summarize, we observe that the negative differential distributions arising from the 
negative Sudakov logarithms yield unphysical singular components present in our hxed-order 
result. When sufficiently broad bins are chosen, these singular components integrate out to 
yield a positive dehnite binned result. In an actual experimental analysis, this requirement 
of broad bins arises naturally given the hnite detector resolution. 

Our use of bins to regularize the differential distributions will have negligible impact 
on the experimental analysis because our bin size is small compared to the experimental 
detector resolution, even more so because the detector acceptance is a smooth function in 
terms of the set of kinematic variables. In particular, while the exact choice of bins is subject 
to any experimental analysis, the geometry of typical neutrino detectors as rehected in the 
acceptance functions of Fig. ^demonstrates the effective experimental bin size in z and rj is 
comfortably large enough for the purpose of regularizing the differential distributions with 
the binning technique. 

To conclude this Section, let us note that in the lowest- 2 ; curves of Fig. (|^), one observes 
~ apart from the forward peak at rj < r^max typical of an hFs' c event - a broader backward 
plateau stemming from the scattering off a collinear (in the boson-target CMS) charm quark. 
For the hxed target kinematics under consideration, this backward peak is seen to be strongly 
suppressed compared to the forward peak. This observation gives us conhdence that one is 
actually measuring the scattering off strange quarks, with only a small dilution from charm 

As an important cross-check on these results, we verify that the binning rjx z = Ixn (rapidity integrated 
out) reproduces the integrated results in the literature |]l 6 |,p^ 

Due to our orientation of the z-axis along the target direction we denote jymax = max (I 77 I) = max (— 77 ) 
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Figure 3: Binned differential distributions in the charm rapidity, r]. Shown are results for 
a hne binning (100x100) in t] x z (left), and a broad binning (right) of 10x5. The hne- 
binned results refer to 2 ;-bins centered around 0.105, 0.605 and 0.955, respectively, while the 
broad-binned results refer to 0.1, 0.7 and 0.9. 


quarks. Hence, this verihes the opposite sign dimuon data provides a direct determination 
of s{x, Q‘^) at NLO level. 


5 Conclusions 

We have presented a fully differential NLO calculation of the neutrino-induced DIS charm 
production process. This calculation is an essential ingredient for a complete analysis of 
the dimuon data, and will allow a precise determination of the strange quark PDF. We 
have demonstrated that by binning the data appropriately, we can interface the theoretical 
calculation (containing ^-functions and “plus-distributions”) directly to the experimental 
Monte Carlo analysis program. We observe the enhancement of the Sudakov logarithms 
at the phase-space boundaries, and verify that these can be controlled with this binning 
method. 

The fully differential distributions obtained here allow charged current neutrino DIS 
experiments to use the complete NLO QCD result in the Monte Carlo data analysis.^ These 
tools will allow us to extract the strange quark PDF from the dimuon data at NLO; this 
^^These results are available as a FORTRAN code, DISCO. 
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Figure 4: Binned distribution in 2 : (left; rj integrated out) and 77 (right: 2 : integrated out), 
each for a hne and a broad binning. 

information should prove crucial to resolving the unusual behavior of the strange quark in 
the proton. 
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